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Bipartite composite boson (quasiboson) systems, which admit realization in terms of deformed 
oscillators, were considered in our previous paper from the viewpoint of entanglement character- 
istics. These characteristics, including entanglement entropy and purity, were expressed through 
the relevant deformation parameter for different quasibosonic states. On the other hand, it is of 
interest to present the entanglement entropy and likewise the purity as function of energy for those 
states. In this work, the corresponding dependencies are found for different states of composite 
bosons realized by deformed oscillators and, for comparison, also for the hydrogen atom viewed as 
composite boson. The obtained results are expressed graphically and their implications discussed. 



I. INTRODUCTION 

Composite bosons or quasibosons or cobosons, as 
non-elementary systems or (quasi-)particlcs built from 
two or more constituent particles, are widely encoun- 
tered [1-7] in diverse branches of modern theoreti- 
cal (quantum) physics. Among quasibosons there are 
mesons, diquarks/tetraquarks, odd-odd or even-even nu- 
clei, positronium, excitons, coopcrons, atoms, etc. In 
the present work we focus on the case of bipartite (two- 
component) composite bosons. Their creation and an- 
nihilation operators can be given through the typical 
ansatz, A\ = $£"0+ 6+ , A a = $0X0^, where 
aj t and b\, are the creation operators for the (distinguish- 
able) constituents, which can be either both fcrmionic or 
both bosonic. In [8, 9] it was shown that the composite 
bosons of particular form (i.e. those that involve appro- 
priate matrices can be realized, in algebraic sense, 
by suitable deformed bosons (deformed oscillators). 

As known, among the measures characterizing degree 
of entanglement or correlation between the entangled 
constituents in a quasiboson there are Schmidt rank, 
Schmidt number, concurrence, purity, and the entangle- 
ment entropy. The latter two are especially important 
in the context of (theoretical and experimental) quan- 
tum information research, quantum communication and 
telcportation [1, 10]. 

It is very important to know how the change of sys- 
tem's energy influences the quantum correlation and/or 
quantum statistics properties of the system under study. 
As it is known, the characteristics of the entangle- 
ment between constituents of quasiboson, which measure 
bosonic quality of quasiboson [11-14], and their energy 
dependence are of importance in quantum information 
research: the quantum communication, entanglement 
production [15], quantum dissociation processes [16], par- 
ticle addition or subtraction in general and in teleporta- 
tion problem [17, 18], etc). Knowledge of the energy 



dependence of witnesses of quantum correlations e.g. en- 
tanglement entropy or purity allows one to relate these 
latter to the energy level of excitation that can be mea- 
sured in the experiments, see e.g. [16, 19, 20]. 

In other words, the energy of a quasiboson differs from 
the energy of the respective ideal boson by a term which 
depends on the quasiboson's entanglement - the measure 
of deviation from bosonic behavior. All this motivates 
to study the energy dependence of the entanglement en- 
tropy and other witnesses of entanglement. In the present 
work we analyze interconnection between the energy of 
system (state) and such main two entanglement charac- 
teristics as entanglement entropy or purity. Note that 
the relationship between the entanglement and energy 
for composite bosons was discussed in [21, 22], for qubits 
in [23, 24], and for spin systems in [25]. 

For those composite bosons realizable by deformed os- 
cillators it is possible, as shown in [26], to link directly 
the relevant parameter of deformation with the entangle- 
ment characteristics of the composite bosons. Namely, 
the characteristics (or measures) of bipartite entangle- 
ment with respect to a- and 6-subsystems, see the above 
ansatz, were found explicitly [26] for single composite bo- 
son, for multi-quasiboson states, and for a coherent state, 
corresponding to the quasibosons system under study. 

Among the above mentioned entanglement characteris- 
tics the entanglement entropy S en t certainly is of primary 
interest. Therefore in this work main attention is de- 
voted to finding the explicit dependence of entanglement 
entropy S en t on the energy E of the quasibosons system 
i.e. of the corresponding state. Present paper further de- 
velops the findings of [26] : we take the composite bosons 
system as being realized in terms of independent-modes 
deformed oscillators with the quadratic 1 structure func- 
tion ip(n) = ^1 + tT^Jn — e|-« 2 , where e = +1/ — 1 for 

fermionic/bosonic constituents respectively. Our analysis 
here is performed for the states considered as the exam- 
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1 As proven in [8, 9] this is the only possibility in case when the 
both constituents are pure fermions (or pure bosons) 
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pies in [26], and also for the hydrogen atom as an inde- 
pendent example. The obtained dependences S en t(E) of 
the entanglement entropy on energy are shown graphi- 
cally for a few values of the deformation parameter /; 
one of the cases is compared with the situation emerging 
for the hydrogen atom. 

Analogous treatment, although in a shorter fashion is 
also performed for the purity-energy dependence. About 
the structure of the paper: in the next Sec. II we sketch 
some facts necessary for what follows; main results on 
the energy dependence of the entanglement entropy of 
(multi-)quasiboson states are presented in Sec. Ill and 
IV, whereas similar treatment for the purity witness of 
bipartite entanglement is given in Section V. The pa- 
per ends with discussion of the obtained results and 
most interesting from our viewpoint physical implications 
(sec. VI). 



II. PRELIMINARIES 



is viewed as bipartite entangled with respect to a- and 
^-subsystems. The degree of entanglement can be mea- 
sured by such well-known characteristics as Schmidt 
rank, Schmidt number, purity, entanglement entropy and 
concurrence, see e.g. [1, 10] for their definition. 

For the entanglement entropy in the case of one com- 
posite boson we obtain [26] , 

2 

S^t = ln(m) = In-, (8) 
whereas for the multi-quasibosonic states (7), see [26], 
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As already mentioned, we deal with composite bosons, 
which are realized by mode-independent system of de- 
formed bosons (deformed oscillators) given for one mode 
by the structure function <p{n). That means that alge- 
braically the quasiboson operators A a , A^ a and the num- 
ber operator N a satisfy on the states the same relations 
as the corresponding deformed oscillator creation, anni- 
hilation and occupation number operators: 



[A a ,A\] = S af3 (^(N a + 1) - (p(N a )), 



[N a ,A 



da/lA^, 



[N a ,A 



-SafiAp, 



(1) 

(2) 
(3) 



where the Kronecker deltas reflect modc-indcpcndencc. 
Such a realization is possible, see [8, 9], only when the 
structure function ip{n) involving discrete deformation 
parameter / is quadratic in n, namely (recall that e = 
±1) 



tp(n)= [l + e- )n-e-n z , f = — , m e IN, (4) 
V 2 / 2 m 

whereas the matrices $ Q , are of the form: 

= C/i(da)diag{0..0,v / 772C/ Q M,0..0}[/ 2 t (4). (5) 

Note that the state of one composite boson 

I^HE^K)*}, K) = at |0>, \K) = fct |0>, (6) 
in/ 

see ansatz, is in general bipartite entangled with respect 
to the states of two constituent fermions (or two bosons); 
likewise, the state describing many composite bosons, 

1*) = ]T *({n 7 })(4 1 )^ • ... • (4J"-|0) (7) 

{n T } 



III. ENERGY DEPENDENCE OF THE 
ENTANGLEMENT ENTROPY 

In order to find the energy dependence of the entan- 
glement entropy we need the expression for the Hamil- 
tonian of the composite boson system. Different choices 
are possible here, but, since quasibosons in our approach 
are realized by means of deformed oscillators, we adopt 
the corresponding Hamiltonian of the same form as e.g. 
in [27, 28]. That is, we take the following Hamiltonian 
of deformed oscillators (deformed bosons) which provide 
realization of the composite bosons: 



H = H \hu a {v{N a ) + <p{N a + 1)). 



(10) 



a. Single composite boson (quasiboson) case. As our 
first example, consider the system which consists of single 
composite boson. For the entanglement entropy in this 
case we have [261 



Seat - In-. 



(11) 



The expression for the energy of one composite boson as 
follows from (10) along with (4), is 

E = ^(1) + p(2)) = - el) = - ±). 

(12) 

Then for the entanglement entropy characterizing single 
composite boson we find 

S ent . = In 



3 E_ 

2 hu 
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The plots corresponding to eq. (13) arc presented 
on Fig.l, Fig. 2. Note the important feature of the 
opposite behavior (increasing vs decreasing) of the en- 
ergy dependence in the case of fermionic constituents 
with respect to the case of bosonic constituents. In the 
both £ = ±1 cases the entropy S en t goes to infinity for 
the energy E = which implies maximal entangle- 

ment between constituents. In this case the constituents 
(fermionic or bosonic) become most tightly bound within 
a quasiboson, and the quasiboson is most close to pure 
boson. On the contrary for E = ^hoj, e = +1, and 
E = §7kj, £ = — 1, the entanglement entropy S cn t = i.e. 
the constituents are unentangled. From physical view- 
point, in this case the constituents are in fact unbound. 



SIF) 




FIG. 1: Dependence of the entanglement entropy S cn t on the en- 
ergy E a for a single composite boson in the case of fermionic com- 
ponents i.e. at e = +1. 
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FIG. 2: Dependence of the entanglement entropy S ont on the en- 
ergy E a for a single composite boson in the case of bosonic com- 
ponents i.e. at e = — 1. 



b. Hydrogen atom as quasiboson. It is of interest to 
consider the hydrogen atom which constitutes a compos- 
ite boson (entangled with respect to proton and electron). 
In this case, however, the relevant matrices are not 
of the form (5), therefore if it was realized (this is an open 
problem) by a deformed boson the latter should be dif- 
ferent from the type mentioned above. So, the creation 
operator for the hydrogen atom with zero total momen- 
tum and quantum number n can be written in second 
quantization formalism (with discrete momenta) as 2 



A 0n — 



V 



(14) 



where ap and 6^L p arc the creation operators for electron 
and proton respectively taken with opposite momenta; V 
is large enough confining volume for the hydrogen atom. 
The momentum-space wavefunction (f> pn is determined 
by the Schrodinger equation: 



^Tr/l) 3 / 2 



e* pr </>„(r)d 3 r; 



h 2 V- 

2m 



4> n (r) + U{v)4> n {r) = E n cj) n (r). 



The expression for the Hydrogen wavefunction in the mo- 
mentum representation is given as [7] 

e**"*> f {2l + l){l-m)\ ^ nm/ ^ l +H\ 

<Ppnlm= ,„ , 1/2 on , -T\ Pi (COS p ' 



(2iry/ 2 V 2{l+m)\ 
'n(n-Z-l)!\i/2 £ 



m ( e-i 

(n+l)\ ) (£2 + 1 )I+2~n-J-l 



rn\n—l-l)\\ 
\ (n + l)l J 



^(fS). a*) 



where P" 1 is the associated Lcgcndrc polynomial, 
Ci~L 1 ;_i(---) is Gcgcnbaucr polynomial, £ = (2n/jh)p, 
7 = Z/na . 

The expansion (14) can be viewed directly as the 
Schmidt decomposition for the state j4p n |0) with Schmidt 
coefficients A p = ^ 2 ^L — 4> pn - Then the entanglement 
entropy for the hydrogen atom is given by the relation 



pi In |A P | 



Sent — — | ^ 
P 



(2nh) 3 
V 



l^pnl 2 



(16) 



where the first equality is nothing but the definition of 
the entanglement entropy [1]. 

Calculation of the expression (16) is moved to ap- 
pendix. Performing derivation we obtain the following 
result 



2 Note that similar ansatz is used for the excitonic creation oper- 
ators, see e.g. [5, 6] 
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-G^lnG^x), (17) 
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where G„,(x) = (1 - x 2 ) l (l - x) 4 (C (x)Y . Let 
consider the simplest case when the quantum numbers 
I = and to = 0. For these values, 



S m t = In 
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s. 



(0) 



Ini. (18) 



Making replacement x = cos a in the integral in (18), 
and using the formula Cl (cos a) = sm (" a ) we infer: 

o n— IV ) sin a ' 

7T 

(0) *i 2 f n 

Sent = Sent ~ In [47m ] / da(l — cos a) sin (na)- 

7T J 



4 sin 2 (na) 



In 



(1 — cos a) 



. 2 
sm a 



(19) 



From the well-known expression for the energy of H- 
atom, E = — Ky/n 2 , we have n = -J—Ky/E. Substi- 
tuting this in (19) we finally obtain 



? (o) 

^cnt 



In 



AS(E) = S ent (E) - S, 
2 f 

/ da(l — cosm sin 

7T J 
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(20) 



The derived energy dependence is shown graphically on 
Fig. 3. As seen, the character of the energy dependence 
here essentially differs from that of the single quasibo- 
son (two-fermion composite) case above, see Fig.l. Main 
reason for distinction lies in that the matrices of 
composite bosons realized by deformed oscillators are dif- 
ferent from the corresponding matrices of hydrogen atom 
given by its wavcfunction <p pn . From the physics view- 
point this implies that the effective interaction between 



the constituents in the above quasiboson is different from 
the Coulomb interaction within hydrogen atom. 

Of course, it would be useful to perform the analysis of 
H-atom system by taking into account the fact that the 
proton, in its turn, also has composite structure (three- 
quark system), thus differing from fundamental (elemen- 
tary) fcrmionic entity. 



E/Ry 



AS(E) -10 



FIG. 3: Dependence of the entanglement entropy AS = S en t- 
from (20) on the energy E a for Hydrogen atom. 
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Let us recall once more that the example of H-atom, 
treated as quasiboson, is included here for comparative 
purpose only. We suppose however that some realization 
(by deformed oscillators) for hydrogen atom as compos- 
ite quasiboson does also exist, though it may be rather 
complicated and non-algebraic. In fact, this task is to be 
solved in two stages: first we have to construct a realiza- 
tion of 3-fcrmion (3-quark) state for the proton subsys- 
tem by some deformed fermionic oscillator, and, second, 
of a quasi-boson formed from the (protonic) quasifermion 
and one more fermion, i.e. electron. We hope to solve 
this very involved problem in near future. 



IV. ENTANGLEMENT ENTROPY VS ENERGY 
FOR MULTI-QUASIBOSON SYSTEM 

Now examine the case of multi-quasiboson states. Tak- 
ing into account the Hamiltonian (10), the total energy 
of the system (at mode-independence) is expressed as 



E 



^2 fkVa 



f 2 

£ 2 n « 



(21) 



a. Qusiboson Fock state. Let us find the entangle- 
ment entropy as function of energy for the normalized 
Fock state of n a quasibosons, [0(n, ct )!] _1 / 2 (Aj l ) Tl ° |0), in 
a fixed mode a. The entropy of entanglement between a- 
and ^-subsystems for the two values of e equals respec- 
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tively, see [26], to 



>cnt|e=+l 



in a 



2//' 



So 



in a 



2/f+n a -V 



(22) 



The latter dependencies S cn t = S en t(n a ) are shown in 
Fig. 4 and Fig 5. 

By inverting eq. (21) we have the dependence of the 
occupation number n a of quasibosons in ath mode on 
the corresponding energy E a of quasibosons: 



K(E a ) 



l±Jl-2ef( 



l) 



2e/ 



Substitution of this expression in (22) leads us to the 
two-branch form of the concerned dependence S^ t (E a ) 
for the case e = +1: 



^t(- B «)|e=+l = ln ( C 2// 



(23) 



where J^j < ^ for ^-branches, > \ for S* cnt - 
branch. For e = — 1 we have single monotonous branch: 

1 [- v /l-/+2/B Q /^„-l]/(2/) 



Sent(S Q )| 6= _i=ln C. 



2/ f-l+[^Jl-f+2fE a /hu a -l]/(2f). 



(24) 



where ^ > ^ . The corresponding functions are presented 
graphically in Fig. 6 and Fig. 7. 
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f=0.05 " 



f=0.2 



f=1.0 



FIG. 4: Dependence of the entanglement entropy Sent, see (22), 
on the number of quasibosons n a for one-mode multi-quasibosonic 
system: the case e = +1 of fermionic components. 

b. The state with one quasiboson per mode. Now let 
us turn to the Example 2 from [26] . In this case the qua- 
sibosons are all in different modes, i.e. the quasibosonic 
system is in the state 

I*) =A\ 1 ■ ... • 4J0), 7*^7j, i¥=3, i,j = l,...,n. 
For the entanglement entropy, for e = ±1, we have 

2 



nln(m) = nln 



(25) 
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f=0.05 ■ 



£=0.2 



f=1.0| 



FIG. 5: Dependence of the entanglement entropy S cn t, see (22), 
on the number of quasibosons n a for one-mode multi-quasibosonic 
system: the case e = —1 of bosonic components. 
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■ f=0.1 



f=0.33 



f=1.0 



FIG. 6: Dependence of the entanglement entropy S en t, see (23), on 
the energy E a for one-mode multi-quasibosonic system: the case 
€ = +1 of fermionic components. 



The energy of the system depends on the dispersion re- 
lation of uJ lj as function of jj. Taking it in linear (in 7j) 
form, namely w 7j . = ujq + (7, — 7i)§^, and also using (21) 
and n 7j = 1 , we arrive at the following expression for the 
energy: 



E =^-^-(hw + ^hujn(n-l)). (26) 



Solving the latter yields n as function of energy, namely 



n(E) 



1 8uj 
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(27) 
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f=0.2 



f=1.0 



FIG. 7: Dependence of the entanglement entropy S en t, see (24), on 
the energy E a for one-mode multi-quasibosonic system: the case 
e = — 1 of bosonic components. 



FIG. 9: Dependence of the entanglement entropy Sent on energy 
E a for multi-quasibosonic system with one quasiboson per mode: 
the case e = — 1 of bosonic constituents. 



where Sui = §^7. Then (25) yields 



Scnt(E) ■ 



1 1 1 5uj 1 / 1 -1 1 Su \ 2 1 a Su 1 E 
~*~ 2 u V \ 2u> ) ' Wo 3-ef too -, 2 

: In — 



Suj/ujq " f 

(28) 

Like in the previous case we obtain the corresponding 
plots which are now placed in Fig. 8 and Fig. 9. 
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f=0.05 f=0.2 



f=1.0 



FIG. 8: Dependence of the entanglement entropy Sent on energy 
E for multi-quasibosonic system with one quasiboson per mode: 
the case e = +1 of fermionic constituents. 

c. Coherent state of quasibosons. As our last exam- 
ple consider the coherent state of composite bosons sys- 
tem in crth mode, see Example 3 in [26]: 



= e-^ 2 /2fi 



1\A\* 



4 TO 



{z/2) m - 1 
, 2 = 2^1^1, 



where I m -i(z) is the modified Bcsscl function of order 
to — 1 . For mean energy of the system in this state we 
have 

E a = (* a \±Hu a [<p(N a ) + ^{N a + l)]|* a ) = ^hw a \C\ 2 - 

00 I A\2n 1 00 I A\2n 
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1 1 n— 

= ^ a (|^ + l/2) + i« Wa i W ^ln^ = 



\A[- 



2 Y TO ^m-l(^) 2 TO 



- + 



Sw a — (to— 1) 



2] 1 ^ , fiWal^j J m {2^i\A\ )+I m -2 (2y/m\A\) 



2m J \frn 



2J m _i(2v^|^|) 



(30) 



The entanglement entropy for the coherent state (29) is 
given, see [26], by the formula (recall that to = 4): 



A n 



\* a ) = C(A;m)Y,^y(Air\0) 



n=0 



(29) 



n=0 



n+m— 1 



(n!) 2 (C» +m -i) 2 



C 2 (|^| 2 TO) f 



(31) 
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Hence we have nothing but the dependence of S en t on E a 
in parametric form (unfortunately, we cannot solve (30) 
for |.4|, here \A\ being the parameter, in order to insert 
the solution into (31); that is why we merely use para- 
metric presentation of the S en t = S cn t(E) dependence). 
The plot of this dependence is given in Fig. 10. 
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E/hm 
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FIG. 10: Dependence of the entanglement entropy Sent in (31) on 
the mean energy E a , see (30), for quasibosonic coherent state. 



OTHER MEASURES (WITNESSES) OF 
ENTANGLEMENT VS ENERGY 



There exist some other widely used witnesses of entan- 
glement: Schmidt number, Schmidt rank, concurrence, 
Schmidt number K or its inverse P — l/K termed pu- 
rity [1, 10]. Energy dependence of these entanglement 
witnesses, Schmidt rank, concurrence and purity have 
somewhat simpler form and can be calculated in a simi- 
lar way using explicit formulas from [26]. 

Since such entanglement characteristics as purity are 
exploited in connection with the issue of entanglement 
creation in scattering processes [15] and others [17, 23], 
let us pay some attention to P. 

For the entangled system consisting of one quasiboson 
the purity in [26] was found to be connected with the 
deformation parameter m = 4 in a simple way: 



k 



(32) 



Then, the energy dependence for purity in case of a 
single composite boson readily follows by combining (32) 
with (12) that gives: 



P 



f 
2 



1/3 
e V2 



tkvJ 



V2 husJ 
fE_ 3 
\~hhj ~ 2 



, 1 E 3 

1, 7^<^<7^, 



e = -1, 



flLO 

3 E 5 

- < — < -. 

2 ~ hui ~ 2 



(33) 



Thus, the dependence of purity on energy is linear for 
both e = +1 and e = — 1. Notice however inverse to each 
other, falling vs raising, behavior of purity with changing 
energy in these cases of fcrmionic vs bosonic constituents. 

In a similar way, for purity in the case of single-mode 
multi- quasibosonic Pock states on the base of [26] we ob- 
tain 

P^(E a)U = +1 = (34) 



Pent(£a)|e=-l = (C 



[y/l-f+2f E a /H^-l]/(2f ) 
2/f-l + ly/l-f+2fE a /hu a -l]/(2fl 



(35) 



the definition intervals being the same as for the entan- 
glement entropy, see (23) and (24). The functions (34) 
and (35) of energy are presented graphically in Fig. 11 
and Fig. 12 correspondingly. Notice the peculiar shape 
of curves in Fig. 11 (non-monotonic behavior, with two 
pieces of monotonicity for each curve). 



P(E) 




■ f=0.1 



f=0.33 



f=1.0 



FIG. 11: Dependence of the purity P on the energy E a for one- 
mode multi-quasibosonic system: the case e = +1 of fermionic 
components. 

Likewise, for multi-quasibosonic states with one quasi- 
boson per mode, using the expression for purity calculated 
in [26] we easily find 



P(E) = ex P (-Sent(£)). 



(36) 



The corresponding plots are now placed in Figs. 13, 14. 

As Figs. 12-14 demonstrate, purity is falling from 
its maximal possible value P = 1 (zero entanglement) 
attained at E < 4^ implying absence of quanta, to 
P = at very large energies of many-particle states. 
A peculiar behavior of purity as function of energy is 
seen in Fig. 11: purity drops from P = 1 with energy 
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| f=0.05 f=0.2 f=l .0 | 

FIG. 12: Dependence of the Purity P on the energy E a for one- 
mode multi-quasibosonic system: the case e = — 1 of bosonic com- 
ponents. 



1 1 




E/fao 

| f=0.05 f=0.2 f=1.0 | 

FIG. 13: Dependence of the purity P on energy E for multi- 
quasibosonic system with one quasiboson per mode: the case e = +1 
of fcrmionic constituents. 



growing to some maximum (the latter is deter- 

mined by the parameter /), then further decreases from 
P f = P(E max (f)) with the energy lowering from £ max (/) 
to smallest values. It is tempting to interpret such exis- 
tence of two regimes as follows: both addition and sub- 
traction [17], [18] of quanta (of quasiboson) can result in 
lowering purity. The two regimes, linked with existence 
of two branches, differ in the starting value of purity P 
(whether it is zero or P(E max )). 



VI. DISCUSSION 

In conclusion, let us make some comments on the ob- 
tained dependencies of entanglement entropy on energy, 




E/to 

| f=0.05 f=0.2 f=l .0 | 

FIG. 14: Dependence of the purity P on energy E a for multi- 
quasibosonic system with one quasiboson per mode: the case e = — 1 
of bosonic constituents. 



and their visualization with the corresponding plots. For 
the state of one composite boson realized by deformed 
oscillator, using the Hamiltonian (10) we find that the 
entanglement entropy monotonically grows with energy 
if the components are fermions, and decreases if the com- 
ponents are bosons (Figures 1 and 2). 

We infer that for larger energies two-fermion quasibo- 
son becomes more tightly bound whereas two-boson qua- 
siboson becomes less bound. In the both cases the energy 
E = 7;hus corresponds to the most entangled quasiboson 
which here shows itself as most close to pure boson. 

If we compare the case of 2-fermion compound with the 
hydrogen atom viewed as two-fermion composite almost 
boson, we observe that in case of iJ-atom the dependence 
of the entanglement entropy on energy shows decreasing 
and thus strongly differs from the two-fermionic quasibo- 
son case (compare Fig. 1 and Fig. 3). The reason may be 
rooted in the specified proton-electron interaction and/or 
in the non-elementary nature of proton, one of the two 
constituent fermions of i7-atom (see also last paragraphs 
in Sec. 3). 

In the case of multi-quasiboson state for a single fixed 
mode and when there are two fermionic components, we 
observe two branches - one decreasing and the other in- 
creasing, see Fig. 6 and the last but one paragraph of this 
section. For the rest of the considered multi-quasibosonic 
states with fermionic components (the fixed one mode 
case, or with one quasiboson per mode, or coherent state) 
the entanglement entropy is monotonously growing with 
energy, see Figs. 7-10, while purity is falling monotoni- 
cally as in Figs. 12-14, or with some peculiarities (two 
regimes or branches of monotonicity, see Fig. 11, end of 
Sec. V, and the last but one paragraph in this last Sec- 
tion) . 

What about the role of deformation parameter /? We 
have quite natural feature: the entanglement entropy is 
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rising with decreasing values of /, i.e. with the approach- 
ing to truly bosonic behavior, either for the Fock states 
at fixed mode or for the coherent states. 

For varying energy, from Figs. 7-10 we find: two- 
fermion quasiboson state, multi-quasibosonic states for 
2-fermion quasibosons within a fixed mode, multi- 
quasibosonic states with one quasiboson per mode, and 
the coherent ones are the more entangled the larger is 
energy. This suggests a possibility to enhance the en- 
tanglement (i.e. its entropy) by increasing energy of the 
(multi-) quasiboson state. The states of two-boson qua- 
sibosons show an opposite behavior as they are less en- 
tangled for larger energies. From the physics viewpoint 
we thus have unambiguous relation between the degree 
(strength) of entanglement and say the energy level of 
the considered multi-quasibosonic states. 

For the states with one quasiboson per mode the de- 
pendence S(E; f) at fixed E is non-monotonic with re- 
spect to deformation parameter /, and in the region of 
small energies of the order 40-60 fku the different curves 
given by differing values of / can intersect (as in Figs. 8 
and 9). It would be interesting to learn more physics 
behind such an intersection, and to confirm that. 

In some cases the dependencies obtained above, e.g. 
those in Fig. 1 and Figs. 7-10, can be viewed in the 
context of entanglement production or enhancement, 
(see [15], [29]) and this provides another possible physi- 
cal implication of our results. As seen, the entanglement 
becomes greater with increasing energy (particle addi- 
tion?) for the listed cases. Unlike those ones, in the 
case of two-boson quasibosons the entanglement creation 
is observed when the energy is decreased. That is, when 
the energy of system is lowered (particle subtraction?), 
the entanglement entropy grows. Possibly, this could be 
checked for some physical examples of the treated sys- 
tems, especially from the particle addition/subtraction 
viewpoint [17], [18]. 

Let us make few more remarks on possible experimen- 
tal verification of the obtained results. That may concern 
the dependencies shown e.g. in Fig. 1 and Fig. 2. As for 
the first case (two-fermionic quasiboson) one may con- 
sider electron-electron or electron-hole composites (exci- 
tons). To test the properties for the two-boson compos- 
ite we could take bi-photons or the H-molecule H 2 for 
the corresponding relevant experiments. Besides, multi- 
quasibosonic dependencies presented in Fig. 8 and Fig. 9 
may also be of practical or physical interest. 

At last, let us note the intriguing appearance of "bifur- 
cations" (or existence of two branches) that are in Fig. 6 
and Fig. 11. Which of the branches is physically realized 
could be an intriguing possible issue for verification. Say, 
2fc-electron, 2fc-photon, fc-exciton systems, etc. could be 
used for dedicated experiments. Besides the quasibosonic 
states studied in this work, non-pure e.g. thermal states 
of quasibosons are also of interest for future analysis. 

We hope to extend the above treatment, the ob- 
tained results, and physical implications, to more com- 
plex quasi-boson (or quasi-fermion) systems in our forth- 



coming works, and also to compare with another real 
physical examples (like the if-atom considered above). 
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Appendix A: Derivation of the entanglement 
entropy for Hydrogen atom 

Transforming the sum in (16) into integral (that im- 
plies very large volume V) and substituting the Hydrogen 
wavefunctions from (15), for 5 en t we obtain 
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